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ABSTRACT 

We have solved numerically the general relativistic induction equations in the interior back¬ 
ground spacetime of a slowly rotating magnetized neutron star. The analytic form of these 
equations was discussed in a recent paper (Rezzolla et al 2001a), where corrections due both 
to the spacetime curvature and to the dragging of reference frames were shown to be present. 
Through a number of calculations we have investigated the evolution of the magnetic field 
with different rates of stellar rotation, different inclination angles between the magnetic mo¬ 
ment and the rotation axis, as well as different values of the electrical conductivity. All of 
these calculations have been performed for a constant temperature relativistic polytropic star 
and make use of a consistent solution of the initial value problem which avoids the use of arti¬ 
ficial analytic functions. Our results show that there exist general relativistic effects introduced 
by the rotation of the spacetime which tend to decrease the decay rate of the magnetic field. 
The rotation-induced corrections are however generally hidden by the high electrical conduc¬ 
tivity of the neutron star matter and when realistic values for the electrical conductivity are 
considered, these corrections become negligible even for the fastest known pulsar. 
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1 INTRODUCTION 

Irrespective of the origin of magnetic fields in neutron stars, whether produced by thermoelectric effects active in a thin layer below the star 
surface when the temperature is much above 10 6 K (see, for instance, Wiebicke & Geppert, 1996), or by a dynamo action during the earliest 
stages of the convective motions (see Thompson & Duncan, 1993), or by post core-collapse accretion of fail-back material after a supernova 
explosion giving rise to a neutron star, a secular decay of the magnetic field is expected as a result of the finite electrical conductivity of the 
stellar matter. The theoretical research in this area is intense, pushed on by the observational evidence that magnetic fields in neutron stars 
are decreasing with increasing spin-down age. There is now a general consensus about the possibility of improving the present knowledge of 
the internal structure of neutron stars by using the constraints from observations of the magnetic field decay. This justifies the effort of taking 
into account all of the possible factors which are supposed to play a role during the decay of the magnetic field. 

Particularly interesting within this context are the general relativistic corrections induced by the presence of a strongly curved back¬ 
ground spacetime. These corrections have been investigated by a number of authors (Ginzburg & Ozernoy 1964, Anderson & Cohen 1970, 
Petterson 1974, Gupta et al 1998, Konno & Kojima 2000) and with a number of different approaches some of which are more rigorous 
(Geppert et al, 2000) than others (Sengupta 1995, 1997). In recent related works, Rezzolla et al (2001a, 2001b) have performed a detailed 
analysis of Maxwell's equations in the external and internal background spacetime of a rotating magnetized conductor. As a result of this 
analysis, it was possible to show that in the case of finite electrical conductivity, general relativistic corrections due both to the spacetime 
curvature and to the dragging of reference frames are present in the induction equations. Moreover, when the stellar rotation is taken into 
account, each component of the magnetic field is governed by its own evolutionary law, thus removing the degeneracy encountered in the 
case of nonrotating spacetimes. The purpose of this paper, which is the natural extension of the work in Rezzolla et al (2001a, hereafter paper 
I), is to quantify the general relativistic effects related to rotation on the evolution of the magnetic field. We have therefore solved numerically 
the general relativistic induction equations derived in Paper I for a relativistic polytropic star with different values of the rotation period and 
of the electrical conductivity. Each of the several calculations performed here benefits from the consistent solution of the initial value problem 
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for a magnetic field which is initially permeating a perfectly conducting relativistic star. This approach avoids the use of artificial initial data 
and provides a more accurate solution of the induction equations. 

Overall, our results show that the rotation of the star and of the background spacetime introduce a decrease in the decay rate of 
the magnetic field. In general, however, the rotation-induced corrections are hidden by the high electrical conductivity of the neutron star 
matter and are effectively negligible even for the fastest known pulsar. Also in the absence of rotation, the spacetime curvature introduces 
modifications to the evolution of the magnetic field when compared with the corresponding evolution in a flat spacetime. These modifications 
depend sensitively on both the metric functions of the interior spacetime and on the radial profile of the electrical conductivity. In the case the 
star is modeled as a polytrope and the electrical conductivity is assumed to be uniform in space and time, the spacetime curvature generally 
increases the decay rate of the magnetic field as compared to the flat spacetime case, with this increase being dependent on the compactness 
of the star. 

The paper is organized as follows: in Section 2 we discuss our treatment of the internal structure of the star in the limit of slow rotation. 
Section 3 is devoted to the solution of the induction equations derived in paper I, with some emphasis on the numerical aspects and in 
particular on the initial value problem. We show our results in Section 4, whereas Section 5 contains the conclusions. Throughout, we use 
a space-like signature (—, +, +, +) and a system of units in which G = c = Mq = 1 (However, for those expressions of astrophysical 
interest, we have written the speed of light explicitly.). Partial spatial derivatives are denoted with a comma. 


2 STELLAR STRUCTURE 


The background metric of a stationary, slowly-rotating star at first order in the angular velocity Cl, is given by 

ds 2 = — e 2<i ^dt 2 + e 2A ^ r) dr 2 — 2 u>(r)r 2 sin 2 ddtdcp + r 2 sin 2 9d<f > 2 , (1) 

where io(r) is the angular velocity of a free-falling inertial frame. For realistic values of the stellar magnetic field (i.e. B = 10 11 — 10 13 
G) we can neglect the contribution of the electromagnetic fields to the background spacetime geometry and determine the internal structure 
of the star and its interior spacetime after solving the following system of ordinary differential equations (henceforth TOV system, from 
Tolmann, 1939; Oppenheimer & Volkoff, 1939) 
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where p(r) is the pressure, e(r) is the energy density and m(r) is the mass enclosed within r. Once an equation of state has been chosen, the 
TOV system can be solved numerically together with the differential equation for the Lense-Thirring angular velocity ui(r) in the internal 
region of the star 
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I is solved from the centre 


where u> = Cl — u>. After selecting a value for the central rest-mass density, the set of differential equations (g 

of the star until the pressure vanishes, thus determining the radius R. For the integration of eq. (|?j), the solution near the centre of the star is 
simplified if we use the analytic power series expansion lo/uj c — 1 + 87r(e c + Pc)r 2 / 5, valid for r —> 0 and where the label “c” refers to a 
quantity at the centre of the star (Miller, 1977). Since in the vacuum region of spacetime external to the star to(r) = 2 J/r 3 , with J being the 
total angular momentum, we can determine the two unknown quantities J and uj c by imposing continuity of the angular velocity and of its 
first derivative at the surface. 

The interior of the star influences the magnetic evolution either macroscopically, by affecting the metric quantities which enter the 
induction equations, or microscopically, through the electrical conductivity a which, in turn, depends on the star’s temperature and chemical 
composition (see Urpin & Konenkov, 1997; Page el al 2000). Our attention is here mainly focussed on assessing the contribution coming 
from rotational effects in general relativity on the decay of the magnetic field Q. As a consequence, we will neglect the thermal and rotational 
evolution of the neutron star and simply consider a constant in time and uniform in space electrical conductivity. This is an approximation, but 
a necessary one to disentangle the many different effects that intervene in the general relativistic evolution of the magnetic field. Furthermore, 
as will be discussed in Section 4, the assumption of a uniform electrical conductivity does not affect the role of a rotating background 
spacetime in the evolution of the magnetic field. 

We model our relativistic stars as polytropes with equation of state 

p = Kp^' N , (4) 


* It should be mentioned that general relativistic corrections can appeal' also in the constitutive relations of the Maxwell equations, such as in the general 
relativistic form of Ohm’s law (Ahmedov 1999). These corrections are usually negligible in the electrodynamics of relativistic stars and will be neglected here. 
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where p, K, N are the rest-mass density, the polytropic constant and the polytropic index, respectively. As “fiducial” model of neutron star 
we consider a polytrope with index N = 1, polytropic constant K = 100, and central rest-mass density p c = 1.28 x 10~ 3 . In this case, 
the radius R and the total mass M obtained through the solution of the TOV system are respectively R = 14.15 Km and M = 1.40 Mg, 
yielding a compactness ratio rj = 0.29. The rotation period usually chosen for this model is P = 1CP 3 s. 


3 EVOLUTION OF THE INTERNAL MAGNETIC FIELD 

As mentioned in the Introduction, the presence of the stellar rotation lifts the degeneracy found in the case of a nonrotating star (Geppert et 
al 2000) and three distinct induction equations regulate the general relativistic evolution of the magnetic field. In this Section we discuss the 
solution of the induction equations for each of the magnetic field components. The main difficulties encountered in the numerical solution 
are related to the definition of a consistent initial value problem and to the complex nature of the partial differential equations when a 
misalignment between the rotation axis and the magnetic dipole moment is present. In the following we discuss the strategies adopted to 
handle these difficulties. 


3.1 The Relativistic Induction Equations 

The induction equations for the magnetic field of a slowly rotating relativistic star with finite electrical conductivity have been derived in paper 
I and we briefly recall them here for completeness. All the measurements are performed in the orthonormal tetrad frame of a “zero angular 
momentum observer” (ZAMO) and we assume that the spatial components of the magnetic field four-vector in this frame are solutions of the 
Maxwell equations in the separable form 


B r {r,9,<t>,X,t) = F{r,t)'Sf 1 {9,<(>,x,t) , 

(5) 

B 9 [r,9, <t>, Xi t) = G(r,t)'5/ 2 (9,(j),X,t) , 

(6) 

B*(r,6,(t>,X,t) = H(r,t)'Sr 3 (<l>,X,t) , 

(V) 

where F, G, H and 4G, '^ 2 , 4^3 account for the radial and angular dependences, respectively. Here, \ is the inclination angle of the stellar 
magnetic dipole moment relative to the rotation axis and the time dependence in F, G, H is here due to the fact that we are not considering 
an infinite electrical conductivity but are allowing the magnetic dipole moment to vary in time. 

At first order in fl, the angular eigenfunctions 4^ are not affected by general relativistic corrections and assume the flat spacetime 
expressions 

4>i = cos x cos 6 + sin x sin 9 cos A(f) , 

(8) 

4-2 = cos x sin 9 — sin \ cos 9 cos A(f) , 

(9) 


q >3 = sin x sin A(t) 


GO) 


where A (t ) = </> — fit is the instantaneous azimuthal position (see Fig. 1 of paper I). Assuming that the contribution of electric currents are 
negligible^] the general relativistic evolution equations for the radial eigenfunctions F(r, t), G(r, f), H (r, t ) are 
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t Because the magnetic field decay is studied on timescales that are much longer than the electromagnetic wave crossing time, this is a very good 
approximation. 
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Together with the evolution equations ( |n| ) - (|IJ), the scalar functions F, G, and H also satisfy the constraint condition of zero-divergence 
for the magnetic field 


j^(r 2 i ? ) r + 2e A rG] sin 9 (cos x cos @ + sin \ sin 6 cos A) + e A r (FI — G) sin x cos A = 0 . (14) 

A rapid look at equations (]ll]) - t [Ti| ) shows that in a rotating spacetime the evolution of the poloidal and toroidal components are correlated 
and that an initially purely poloidal magnetic filed can gain a toroidal component during its evolution and vice-versa. In the case of a 
nonrotating star, on the other hand, the three induction equations f|TT|) - © are not independent and the magnetic field evolution is described 
by a single scalar function: F (see Geppert et al. 2000). 


3.2 Strategy of the Numerical Solution 

The numerical solution of equations (JTTJ) — (JT^) is simplified if done in terms of the new quantities 

F = r 2 F , 


G = e rG . 


(15) 

(16) 


77 = e*rH , (17) 

which, when the inclination angle \ is nonzero and the electrical conductivity is uniform, allow us to rewrite eqs. ( |Tl| ) - (|T^) schematically 
as 
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Explicit expressions for the set of coefficients fi,gt, hi can be found in Appendix A. For \ 7 ^ 0, the coefficients fi, gt, hi have terms which 
are time-dependent trigonometric functions of Fit and, as a result, each of the eqs. (jTi[) - (|2c|) is not a simple parabolic equation describing 
a pure diffusive phenomenon. In addition to a secular Ohmic decay, in fact, there will be a periodic modulation produced by the rotation of 
the star. This is evident if we look, for instance, at the coefficient /1 in Appendix A and which is given by the sum of two terms. The first 
one is the constant “diffusion” coefficient responsible for the decay on a secular timescale. The second term, on the other hand, represents 
the correction due to the stellar rotation. The periodic modulation is produced by the trigonometric function tan A and varies therefore on the 
dynamical timescale set by the rotation period of the star, P. The presence of these periodic terms spoils the parabolic character and makes 
the set of eqs. (|T^) - (^(j) a mixed hyperbolic-parabolic one. 

Although the integration of eqs. (jF 8 ]) - ( |2ti| ) is complicated in the general case, we are here favoured by the fact that all of the terms 
proportional to fl or to ui (i.e. all of the terms directly related to the stellar rotation) scale like a ~ 2 and that the electrical conductivity in 
realistic neutron stars is very high, ranging in the interval 10 21 — 10 28 s -1 . As a result, the star’s rotation period is about twenty orders of 
magnitude smaller than the secular timescale and can be ignored in the numerical solution of the equations. In practice then, we set all of the 
periodic time-varying terms to be constant coefficients and solve the set of eqs. (|l 8 j) - ( p(i[ ) as a purely parabolic system. In this way we can 
capture the secular decay without having to pay attention to the high frequency modulation. In Section 4, where we discuss the results of the 
numerical integration of the induction equations we will also comment on the validation of this procedure. 

Another important aspect of the numerical solution is the use of the zero-divergence constrain equation (^j). We do not need, in fact, 
to integrate in time all of the eqs. but can restrict the evolution to two of them and obtain, at each timestep, the remaining 

unknown radial eigenfunction from the solution of the constraint equation ©• Adopting this strategy in the numerical solution reduces the 
computational costs and, most importantly, enforces a constrained solution at each timestep. 
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Having three induction equations, we can follow the decay of each component of the magnetic field separately. The physically relevant 
quantity is however the modulus of the magnetic field which, in the locally flat spacetime of the ZAMO observer, is simply given by 
|J3| = [( B r ) 2 + (H®) 2 + ( B < ^) 2 ] 1 / 2 . The evolution of this quantity, evaluated at the surface of the star, is the one that will be discussed in 
the remainder of the paper. 


3.3 The Initial Value Problem 

The consistent solution of the initial value problem for the general relativistic decay of the magnetic field in a rotating neutron star suffers 
from two difficult aspects. The first one is that at present the initial topology and location of the magnetic field in neutron stars can be only 
argued on the basis of some assumptions, so that the magnetic field can either permeate the entire star, or be confined in a layer close to 
the stellar surface. The first configuration is more plausible if the magnetic field is the final product of a dynamo action amplification (see 
Thompson & Duncan, 1993), while the second field configuration is more realistic in a scenario in which the magnetic field is originated by 
thermoelectric effects (Urpin et al, 1986; Wiebicke & Geppert, 1996). We here focus our attention mostly on the case of a magnetic field 
permeating the entire star, but in Section 4 we also show how the decay of the magnetic field depends on the depth of penetration inside the 
star, when simplified assumptions on the microphysics at the crust-core boundary are made. 

The second difficult aspect of the initial value problem concerns the definition of an initial configuration which is also solution of the 
general relativistic Maxwell equations. A possible approach to this problem is the one proposed by Geppert et al (2000) (but see also Sang 
and Chanmugan, 1987), who have considered the initial magnetic field to be described by Stoke functions that represent, in flat spacetime, a 
class of exact solutions of the induction equation. In this case, the radial eigenfunction F(r) at the initial time can be obtained from 

F(r, t) = B 0 

for t = 0, where Bo is the initial surface magnetic field at the magnetic pole. Because eq. (^Tj) is not a solution of the general relativistic 
Maxwell equations, one expects an initial error being introduced in the solution of the induction equations, but also that this error should 
disappear rapidly as the solution tends to the one satisfying the Maxwell equations. 

To circumvent the problem of an inaccurate solution during the initial stages of the evolution and in order to calculate an initial magnetic 
field which is solution of the relativistic Maxwell equations, we here treat the initial magnetic field as the one permeating a perfectly 
conducting medium. In this case, Rezzolla et al. (2001a, 2001b) have shown that consistent radial eigenfunctions can be obtained after 
solving the following set of equations [see (71)-(73) of paper I] 

F, r + 2e A_ *G = 0 , (22) 
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In particular, combining eqs. 


and ([23]), we obtain a second-order differential equation for the unknown radial eigenfunction F 


d?F d dF „ 2A F „ 

~dr T+ dr^ ~ A ^7i7~ e 7? ~°- 


dr 


(25) 


Equation (|0 can be solved as a two-point boundary value problem after specifying values for the magnetic field at the edges of the 
numerical grid. More specifically, the initial magnetic field at the inner edge of the grid is chosen to be zero both when the magnetic field 
permeates the whole star and when it is confined to a crustal layer. On the other hand, the initial magnetic field at the outer edge of the grid is 
chosen to match a typical surface magnetic field for a neutron star and is therefore set to be Bo = 10 12 G. Once the initial profile for F has 
been calculated through eq. (|25|), the corresponding initial values for G and Ft follow immediately from eqs. ( |22| ) and (0) . As a comparison, 
we have also solved the induction eqs. ([i~s[) (]2o|) using as initial condition eq. (0) and the corresponding eigenfunctions G and Ft again as 
computed from the conditions ( |22| ) and (|24|). 

Fig. [I] shows the initial values for the two different prescriptions and, in particular, with a solid line the initial profile as obtained through 
the solution of the Maxwell eqs. ( [25| ) and with a dashed line the Stoke profile given by expression (fTl]). The noticeable differences between 
the two initial profiles provide a simple explanation of why the use of Stake’s function produces an initially inaccurate evolution (cf. Fig.|^). 

The use of the strategy discussed above for the calculation of the initial value problem clearly requires the solution of an additional set 
of equations but it has the advantage of removing the adjustment of the solution during the initial stages of the decay and provides a more 
accurate estimate of the magnetic field decay. A discussion of this as well as a comparison with evolutions performed with the Stoke function 
will be discussed in Section 4. Finally, another aspect worth stressing is that by using eq. ( ]25| ) we also automatically satisfy appropriate 
boundary conditions at the surface of the star. 
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Figure 1. Possible initial values for the radial eigenfunction F of the radial component of the magnetic field normalized to the value at the surface, shown 
as a function of the radial position in the star. The solid line represents the radial eigenfunction F as obtained from the integration of the Maxwell eqs. (jzij), 
while the dashed line represents the Stoke profile 0 - 


3.4 Boundary Conditions 

In order to correctly solve the induction equations ([Ti]) - it is essential that appropriate boundary conditions are specified both at the 
inner edge of the computational domain as well as at the stellar surface. 

As for the initial value problem, the inner boundary condition imposed during the evolution is that of a zero magnetic field and is applied 
both when the magnetic field permeates the whole star and when it is confined to the crust. In the first case, this choice guarantees a regular 
behaviour of the radial eigenfuctions at the origin, while it reflects the absence of magnetic field below the crust in the second case. The 
evolution of the magnetic field has shown to be quite sensitive to the boundary conditions imposed at the stellar surface, but proper boundary 
conditions can be derived if we assume that there are no electrical currents on the surface and impose a matching between the external and 
the internal solutions of the magnetic field. The radial eigenfunctions F(r), G(r), and H(r) outside the slowly rotating relativistic star have 
been derived in paper I [see eqs. (90) - (92) therein] and are given by 
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4 M 3 
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In N 2 H- 
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H(r) = G(r) , (28) 

where N(r) = (1 — 2M/r )V 2 = e* and /.i is the magnetic dipole moment. Since the constraint expressed by eq. ( p2| ) holds also on the 
stellar surface, we then have 


F r 


+ 2e A_s G(J?) = 0 . 


Moreover, when electrical surface currents are not present, we can use eq. ( |26| ) and ( |27| ) to express G{R) as 

N 2 M\ Rhi N 2 /M + A/N 2 + 1 


G{R) = - — 5 - —=-- '—!■ --- F{R) , 
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where N = N(r = R). Straightforward calculations allow to conclude that 
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Figure 2. Difference in the decay of the magnetic field when a consistent initial magnetic field is used (solid line) or when a Stoke function is used as initial 
condition (dashed line). The inset shows a magnification of the evolution during the first 10 9 yr. Here a = 10 25 s —1 and P = 10 —3 s. See the main text for a 
complete discussion. 


where 11 ( 77 ) is a constant given by 


nfa) 


4 In (1 - nj) + 277(2 - ?;)/(! - r?) 
2 In (1 — rj) + 2rj + r/ 2 


(32) 


with ri = 2 M/R being the compactness of the star. The corresponding boundary conditions for G and H are then easily obtained by means 
of (|§) and @). 

Note that eq. ( [Tij ) coincides with the boundary condition used by Geppert et al (2000) in the case of a static, spherically symmetric 
background geometry. This is due to the fact that, as discussed in paper I, the corrections to the components of the magnetic field enter at 
orders higher than the first one in Cl. Details on the numerical implementation of the surface boundary conditions are presented in Appendix 
B. 


4 NUMERICAL RESULTS 

In order to integrate the set of induction eqs. ( [l8[ )- we have built a numerical code which implements the Crank-Nicholson implicit 
evolution scheme and which provides second order accuracy both in space and in time (see Morton & Mayers, 1994). The accuracy of the 
code has been checked by computing the time evolution of eq. 0 which provides, in a flat spacetime, an exact solution of the induction 
equation. The results obtained indicate that the relative error between the numerical and the analytic solutions over a timescale of three 
Newtonian Ohmic times r 0 hm = 'inlEa/c 2 , is always below 0.5% for the level of grid resolution usually implemented in our calculations. 

Established the consistency and accuracy of the code, we have proceeded to solve the general relativistic induction equations for our 
relativistic rotating star. As mentioned in Section 3.2, if the inclination angle between the rotation axis and the dipolar magnetic moment is 
nonzero, the secular decay has a periodic modulation due to the stellar rotation. We have also discussed that because the decay timescale and 
the rotation period timescale differ for about twenty orders of magnitude, we can neglect the time dependence (which is oc sin A) contained 
in each of the coefficients fi.gi, hi and set the periodic terms equal to an arbitrary constant value. To validate this procedure and verify that 
the periodic modulation does not affect the secular evolution, we have solved the induction equations using different constant coefficients and 
found that the secular results are indeed insensitive to the value chosen for the constant coefficients. We have also followed the solution of the 
complete set of eqs. (i.e. not considering the time-periodic terms as constant) on a timescale which is longer than the dynamical 

timescale but still much smaller then the secular one. Also in this case we have verified that the modulated evolution, which is superimposed 
on the secular one, shows the small decrease corresponding to the secular decay. 

Our discussion of the results starts by comparing the evolution of eqs. ([Hi|) - (^) for the two different prescriptions of the initial value 
problem discussed in Sect. T2 (cf. Fig. |l|) for our fiducial neutron star. Before presenting the results of the comparison, it is useful to discuss 
briefly the subtleties related to the measure of the magnetic field time decay; as will become apparent later, this is an important issue which 
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Figure 3. Relative difference in the evolution of a magnetic field in a nonrotating star Bnonrot and in a rapidly rotating one B ro t. The electrical conductivity 
is here set to be a = 10* s — , while the star been set to have a period P = I () 3 s. 


might lead to seemingly conflicting results. The gauge freedom inherent in the theory of General Relativity allows for the choice of arbitrary 
observers with respect to which the measure of physically relevant quantities is made. The choice of a certain class of observers might rely 
on the mathematical advantages that this class may have, but not all observers are physically suitable observers. Locally inertial observers 
are certainly preferable and in a rotating spacetime, as the one considered here, ZAMO observers represent a natural choice. Of course, there 
is is an infinite number of such observers, each one performing his own measure of the magnetic field decay, so that one should then select 
a specific set of inertial observers on the basis of physical considerations. The results presented in this paper will be referred to a ZAMO 
observer on the surface of the star and at a latitude 9 = n /2. The values of the magnetic field measured by this observer and its time evolution 
can then be converted to the equivalent ones measured by other ZAMOs at different radial and polar positions through simple transformations 
involving the difference in the red-shifts and latitudes. Once the choice of a suitable class of inertial observers is made, it is also important 
that the results of the general relativistic magnetic field decay are expressed using appropriate units. In their work, Geppert et al (2000) have 
quantified the decay of magnetic field in a relativistic constant density, nonrotating star in terms of the Newtonian Ohmic time. As we shall 
show below, while this choice is acceptable for a constant density star, it could be misleading in general. 

The two solutions of eqs. (p~8|) - (jic|) are presented in Fig. ^ and show the decay of the magnetic field , rescaled on a timescale tg = 10 9 
yr. It is interesting to note that while the asymptotic decay rates of the magnetic field are almost the same for the two approaches, a final 
difference emerge. This is because when using Stoke’s function as initial data, the evolution does not satisfy Maxwell’s equations during 
the initial stages (see the small inset in Fig. |2j), but settles onto a constrained solution only after that time. Moreover, the outer boundary 
conditions expressed by eqs. (jicj) and ( pTj ) cannot be satisfied exactly by Stoke’s function and this introduces an additional error. As a result, 
after a time t ~ tg yr, the two solutions differ of about 45%, but this difference does not grow further in time. 

Next, we discuss the effects introduced by the rotation of the star and of the spacetime. In this case it is worth distinguishing the interest 
in finding a general relativistic correction, from the impact that these corrections actually have on the magnetic field decay in a realistic 
rotating neutron star. As discussed in Section 3.2, in fact, the high value of the electrical conductivity in realistic neutron stars tends to make 
the general relativistic corrections due to rotation rather minute. In particular, we have found that when considering an electrical conductivity 
a = 10 25 s' 1 in a rapidly rotating neutron star with one millisecond rotation period, the relative difference in the magnetic field after 15 tg 
yr is only one part in 10 12 . Nevertheless, general relativistic, rotation-induced corrections have an interest of their own and these corrections 
can be more easily appreciated if smaller (and therefore less realistic) values of the electrical conductivity are considered. In Fig. || we show 
the relative difference in the evolution of a magnetic field in a nonrotating star, B nonro t, and in a rapidly rotating one with a millisecond 
period, B ro t. In this case and just for illustrative purposes, an electrical conductivity a = 10 s s _1 has been considered. As can be appreciated 
from the figure, the corrections due to the rotation decrease the rate of decay of the magnetic field and after a few rotation periods, the fastly 
rotating star will maintain a magnetic field which is about a factor of two larger than the one calculated for the nonrotating star. Overall, 
the results obtained indicate that General Relativity does introduce, through the rotation of the spacetime, new corrections to the evolution 
of the magnetic field, slightly decreasing its decay rate. This effect, however, is usually hidden by the high electrical conductivity of the 
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Figure 4. Decay of the surface magnetic field as measured by a ZAMO observer on the surface of the star at a latitude 9 = n/2, expressed on a timescale 
tg = 10 9 yr. The left panel refers to a constant density stellar model and shows an asymptotic decay rate of the magnetic field which is decreasing for 
increasing values of the stellar compactness. The inset in the left panel focuses on the initial stages of the evolution when the decay is larger. The right panel, 
on the other hand, refers to an N = 1 polytropic stellar model and shows an asymptotic decay rate which is increasing for increasing values of the stellar 
compactness. Here the central density is the free parameter determining the stellar compactness. The small inset in the right panel of the figure shows how the 
use of an Ohmic timescale as normalizing unit can lead to erroneous interpretations. 


stellar medium and can be neglected in general. The results discussed above depend also on the inclination between the rotation axis and the 
magnetic dipole moment, with the decrease rate being larger for larger inclination angles. In particular, for x = 7r/2, the residual magnetic 
field after 10 tg yr is smaller of a factor two as compared to the corresponding magnetic field for an inclination x = 0. 

Next, we compare the results of our calculations for a polytropic relativistic star with those for a constant density star. This will provide 
a first qualitative estimate of the importance of the metric functions in the actual evolution of the magnetic field. The results are presented in 
Fig. ^ with the left panel referring to a constant density model and the right one to our fiducial polytropic model. 

In the case of a constant-density star, we confirm the results obtained by Geppert et al (2000) and find that the evolution of the magnetic 
field approaches an exponentially decaying behaviour, with an asymptotic decay rate which is generally decreasing with increasing stellar 
compactness. The inset in the left panel of Fig. ^ shows in more detail the initial stages of the magnetic field decay and allows to appreciate 
that the magnetic field evolution is initially following an exponential decay with decay rates which are quite large but that then reach an 
asymptotic value after about 10 8 yr (Geppert et at 2000). 

In the case of a polytropic star, on the other hand, the results in the right panel of Fig. |dj show a behaviour which is the opposite to 
the one encountered for a constant density model and that the asymptotic decay rate of the magnetic field is increasing with increasing 
stellar compactness. When a uniform electrical conductivity is used, the explanation behind the two distinct behaviours has to be found in 
the deviations that emerge in the internal spacetime for the two stellar models and in particular in the first radial derivatives of the metric 
functions <1? and A [cf. eqs. These deviations produce sensible quantitative differences in the coefficients of eqs. (see the 

Appendix for the explicit form of the coefficients) which are then responsible for the increase in the decay rate. It should also be remarked 
that the behaviour shown in the left panel could be easily reproduced, also in the case of a polytropic model, by means of a suitably defined 
electrical conductivity. In other words, the results presented in Fig. ^ underline that a definitive conclusion on the general relativistic evolution 
of the magnetic field cannot be drawn until a more realistic treatment of the electrical conductivity and of the equation of state is made. 

The inset in the left panel of Fig. ^ can be used to explain the comment made above on the use of relevant normalization units. In the 
inset, in fact, we have plotted the same evolution shown in the main panel but with the time being normalized in terms of the Newtonian 
Ohmic time. Note that when we do so, the overall behaviour is inverted and the decay rate of the magnetic field in now decreasing for 
increasing stellar compactness. This is clearly incorrect and the misleading behaviour is due to the fact that the concept of an Ohmic time is 
a purely Newtonian one and is therefore justified only in a Newtonian context. A more suitable normalizing unit for a nonrotating relativistic 
star would be the general relativistic analogue of the Newtonian Ohmic time: r 0 hm = 47r_R 2 e 2A_# cr/c 2 as can be derived from eqs. (|T^) - 
( ^oj ) in the limit 0 = 0. Using this normalization, we would recover the correct behaviour, with a magnetic field asymptotic decay rate 
generally increasing with stellar compactness. Unfortunately the validity of r 0 h m is limited to nonrotating stellar models only. Because of the 
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Figure 5. Decay of the surface magnetic field when the magnetic field does not penetrate the whole star. The different curves refer to different values of the 
parameter q = R m / R, with f? IN being the inner radius of the stellar shell where the magnetic field is confined. 


difficulties of defining an Ohmic timescale for the induction equations of a relativistic rotating star, we measure the magnetic field evolution 
simply in terms of the time measured by our ZAMO observer. 

Finally, we discuss the differences introduced in the decay of the magnetic field when the latter is confined to a spherical shell between 
an inner radius R m and the surface of the star. In this case, the initial values for the radial eigenfunctions are calculated self-consistently 
along the procedure discussed in Section 3.3. In Fig. |5] we show the evolution of the magnetic field in our fiducial neutron star for different 
values of the parameter q = R m /R. Note that decreasing the volume in which the magnetic field is confined has the effect of increasing the 
decay rate of the magnetic field so that if the initial magnetic field permeates about 90% of the stellar volume (q = 0.5), the residual surface 
magnetic field after 10 ig yr is about a factor thirty smaller than in the case the magnetic field permeates the whole star (q = 0). Although 
our analysis does not take into account the microphysics of the stellar interior and in particular the role played by the chemical composition 
and by the temperature, it confirms the Newtonian results of Urpin and & Konenkov (1997) and those of Page et al (2000), who have shown 
that the magnetic field decay is slower for deeper magnetic field penetration. Because this behaviour mimics the increase in the decay rate 
produced by an increasing compactness of the stellar model, it is essential to be able to determine, prior to observations, the geometry and 
location of the magnetic field within the neutron star and to distinguish the different contributions to the overall magnetic field decay. 


5 CONCLUSIONS 

In a recent paper, Rezzolla et al (2001a) have considered the general relativistic description of the electromagnetic fields of a slowly rotating, 
magnetized and misaligned neutron star. If the stellar medium has a finite electrical conductivity it was shown that the stellar rotation removes 
the degeneracy in the evolution equations for the magnetic field and that three distinct induction equations need to be solved to account for 
the decay of the stellar magnetic field. In this paper we have solved numerically the general relativistic induction equations derived in paper 
I, investigating the effects of different rotation rates, different inclination angles between the magnetic moment and the rotation axis, as well 
as different values of the electrical conductivity. The aim of these numerical calculations is that of quantifying the corrections induced by 
general relativistic effects (both due to spacetime curvature and to the stellar rotation) on the evolution of the magnetic field of a slowly 
rotating neutron star. 

In order to single out purely general relativistic effects from those due to the microphysics of the Ohmic dissipation, we have considered 
a simplified physical description of the neutron star. In particular, the star has been modelled as a polytrope rotating with a fiducial period of 
one millisecond, the electrical conductivity has been considered to be uniform inside the star and we have not included a treatment to consider 
the evolution of the stellar rotation and temperature (see Page et al 2000). On the other hand, special attention has been paid to a consistent 
solution of the initial value problem and we have considered as initial magnetic field the stationary solution of the general relativistic Maxwell 
equations. In this way we have avoided the use of initial magnetic field configurations that are only approximate solutions of the Maxwell 
equations (i.e. solutions of the Maxwell equations only in the limit of flat spacetime). Besides eliminating an initial error during the initial 
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stages of the magnetic field decay, our prescription for the initial value problem also provides a more accurate solution of the Maxwell 
equations. 

The results of our computations have shown that there exist general relativistic, rotation-induced corrections to the evolution of the 
magnetic field. These effects generally produce a decrease in the rate of magnetic field decay. However, their contribution is masked by the 
high value of the electrical conductivity in realistic neutron stars and can be neglected in general. Our calculations also indicate that general 
relativistic effects not induced by the stellar rotation can modify the time evolution of the magnetic field in a magnetized star. Such effects 
are closely related to the properties of the spacetime internal to the star and for a polytropic stellar model with uniform electrical conductivity 
these effects generally increase the decay rate of the field. The validity of this conclusion is however limited. Density gradients are in fact 
expected in a realistic star and these will affect the behaviour of the electrical conductivity which, in turn, will influence the decay of the 
magnetic field. 

Our conclusions are that the general relativistic evolution of the magnetic field in rotating neutron stars can be studied with confidence 
already in a nonrotating background spacetime. However, the role of a curved background spacetime on the decay of the magnetic field can 
be fully assessed only when the details of both a realistic equation of state and of a realistic electrical conductivity are carefully taken into 
account. This will be the subject of future work. 
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APPENDIX A: THE NUMERICAL SOLUTION OF THE INDUCTION EQUATIONS 

In this Appendix we provide the explicit expressions for the coefficients fi, gi , hi appearing in the new form of the induction equations 
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APPENDIX B: IMPLEMENTING SURFACE BOUNDARY CONDITIONS 

This Appendix shows how the surface boundary conditions expressed by equations <0 and (|T|) can be implemented in a numerical code. 
By adopting the standard finite-difference notation in which u" = u(xj, t n ) and assuming a uniform radial grid with J gridpoints, the finite 
difference form of eq. t [TT[ ) is given by 

Fj+i - F?-! = 2U( V )AxF?/R, (Bl) 
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where Ax = Xj — Xj-i, At = t n+1 — t n . The unknown value of Fj +1 , comes after introduction of © into the Crank-Nicholson scheme, 
centered at J; lengthy but straightforward calculation give 
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where a = At/Ax 2 . The are still two unknowns entering (B2), i.e. Gj +1 and Hj +± . However, they represent the external solution and by 


using eqs. I 
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The updated values of G and H now follow immediately from (|3C|) with time evolved value of FJ +1 given by 
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